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1. Introduction 

Among the most important quantities playing a role in the theoretical interpretation of cur¬ 
rent heavy ion collision experiments at RHIC (Brookhaven National Laboratory, USA) and LHC 
(CERN, Switzerland) are so-called transport coefficients: shear and bulk viscosity’s as well as 
heavy and light quark diffusion coefficients. Because of strong interactions, these quantities need 
to be determined by non-perturbative numerical lattice Monte Carlo simulations. This task is a hard 
one, given that numerical simulations are carried out in Euclidean signature, whereas transport co¬ 
efficients are Minkowskian quantities, necessitating an analytic continuation [1]. Nevertheless, the 
problem is solvable in principle [2], provided that lattice simulations reach a continuum limit and 
that short-distance singularities can be subtracted [3]. We present selected results of transport co¬ 
efficient calculations and scale setting in pure SU{3) gauge theory. Scale setting is involved in in 
the proper continuum extrapolation of lattice correlation functions. 

Most of the results concerning scales and diffusion in the SU (3) plasma have been published 
in [4, 5] or will be presented elsewhere [6, 7] . Here we specifically supply additional insight in the 
determination of Tc, which at finite temperature is the most fundamental physical observable. We 
leave out our groups high precision studies of the technical lattice scales ro/a, the Sommer scale 
[8] as well as the Wilson flow scale ^/toja [9]. We will however make use of the data in [10] and 
determine the "amplitude product" Tcy/ff using Tc data and the Euscher et. al. to data set. 


2. Colour electric Correlation Function 


Heavy quarks carry a colour charge and, whenever there are gauge fields present, are therefore 
subject to a coloured Eorentz force. Eike with other transport coefficients the corresponding “low- 
energy constants” are easiest to define af vanishing three-momentum; then the Eorentz-force is 
proportional to the electric field strength. This leads to a “colour-electric correlator” [11, 12], 


Ge{^) 



ReTr 


T)gEKuO)U(T,0)gE,(0,0) 


(ReTr[U(d,0)]) 


(2.1) 


where gEi denotes the colour-electric field, T the temperature, and f7(T2j'^i) a Wilson line in Eu¬ 
clidean time direction. If the corresponding spectral function, pg, can be extracted [1], then the 
“momentum diffusion coefficient”, often denoted by K, can be obtained from 


K = lim 

(B^O 


2Tp^((o) 

(0 


( 2 . 2 ) 


According to non-relativistic linear response relations (valid for M ^ tiT, where M stands for a 
heavy quark pole mass) the corresponding “diffusion coefficient” is given by D = 2T'^ jK. We 
determine stochastic estimates for the function Ge(t) from large scale Monte Carlo simulations in 
finite temperature T « \A5Tc pure gauge SU (3) lattice gauge theory with standard Wilson action 
and for box sizes, that substantially exceed those of earlier simulation results [13]. The largest 
volume is Umax = 48 X 192^ which is simulated with sufficient statistics at j3 = 7.793. Table 1 shows 
run parameters with Astat labeling the number of statistically independent configurations, while r^T 
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Nr 

Ns 

i^conf 

vqT 

6.872 

16 

64 

100 

1.111 

7.192 

24 

96 

160 

1.077 

7.544 

36 

144 

563 

1.068 

7.793 

48 

192 

223 

1.055 


Table 1: A lattice discretized version of the correlator eq. 2.1 is calculated at the given j3 values and box 
sizes for the finite high temperature T « 1.45rc. 



Figure 1: Correlation function ratio GimpiTT)/GnovmijT), see text, with perturbative NLO [19] and NLLO 
calculations. The black band denotes the desired continuum extrapolation result. 


denotes an estimate of the Sommer scale in units of the temporal box extent. We perform quenched 
QCD calculations for the discrete version of the correlation function G£(t). Due to the bad signal 
to noise ratio of the purely gluonic operator and the rapid decrease of the correlation function with 
T we have to use noise reduction techniques i.e., multi-level updates [14, 15] and link-integration 
“PPR” [16, 17]. Furthermore we use a tree-level improvement [18] to reduce cut-off effects and 
a NLO perturbative renormalization factor Zpeni^)- In Fig. 1 the lattice results for Gi,„p{'^T) are 
shown, normalized to 


Onorm('fF) = 


cos^ {titT) 1 

sin^(7rTr) 3sin^(7rTr) 


(2.3) 


Although cut-off effects are visible at small separations and the results become more noisy at large 
distances on the finer lattices, the results on the four lattices allow for a controlled continuum 
extrapolation down to distances around zT ~ 0.05. To obtain the continuum estimate of the corre¬ 
lation function, we perform b-spline or polynomial interpolations for each lattice and at fixed zT 
extrapolate the correlator in 1 /N^. To allow for needed contributions to the spectral function PE(ft)) 
we use the Ansatz 


Pmodel{Oi) = m&X^ApNNLo{0^) +B(0^, - \ 


2T i ’ 

containing three parameters A, B and K, and fit the Euclidean correlator with the form 


/-“dm cosh(I-Tr) 

Gmodel{r)-j^ ^Pmodel{(0) 


(2.4) 


(2.5) 


3 


















Continuum Study of Heavy Quark Diffusion 


T. Neuhaus 


Nx 

Pc 

4 

5.69275(28) 

6 

5.89425(29) 

8 

6.06239(38) 

10 

6.20872(47) 

12 

6.33514(45) 

14 

6.4473(18) 

16 

6.5457(40) 

18 

6.6331(20) 

20 

6.7132(26) 

22 

6.7986(65) 


Table 2: Numerical fc values for pure gauge SU{3) theory with Wilson action. 
2.2 
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Figure 2: Given the critical time extent tcja = Nx as calculated from Table 2 we display the quantity 
ln(fc/a) — 47r^(j3 — 6)/33 as a function of j3. The curve for j3 > 5.8 {Nx > 4) displays the fitted rational 
Ansatz eq. 3.1, see the text. 


in the xT range [0.1:0.5] to the continuum. Under the given procedure we find the momentum 
diffusion coefficient: 


K-/r^= 2.5(4) r«1.45rc. (2.6) 

For a more detailed presentation we refer here to [5]. We will have to investigate the stability of 
the numeric value under various model assumptions for the spectral density. 

3. Deconfinement Position at large temporal Extent 

Using alternate finite size scaling arguments [20] from statistical physics, based on phase 
weight ratios at first order phase transitions we determine the finite temperature transition couplings 
l5c{Nx) for a maximal value Nx = 22. The numerical findings are contained in Table 2. The 
Pc data of Table 2 correspond directly to the circles in Fig. 2 where the temporal time extent 
tcja = Nx in units of a at the SIJ (3) Deconfinement phase transition is suitably logarithmized and 
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Figure 3: For the Nt = 22 and Ns = 64 box we display as a function of j3 the step function s(j3), whose zero 
determines the critical point at = 22 as in Table 2. The zero is marked by a solid circle and horizontal 
error bar and is obtained by the straight line fit, the line in the figure. 


22 X 6r 



subtracted: \n{tc/a) — 4n^{[5 — 6 )/33, and then is displayed as a function of j3. Error propagation 
via 5tc ~ (471^/33) was used. It has become customary, as in case for the Sommer scale ro, 

to interpolate logarithmized QCD scales \n{s/a) with a rational Ansatz [21] 





1 +cip 
1+C3j8 


‘+C2i8-2 

'+C4i8-2^ 


(3.1) 


at values bQ = \\/= 102 /( 47 r)^. We find the fit parameter values 


Cl = -9.5347(2060) C 2 = 22.252(1.053) C 3 = -6.5581(2958) C 4 = 6.547(1.241) (3.2) 

at a reduced Zd o f ~ ’^^^e fit if all data points with Nz >6 are included. See the curve in 

Fig. 2. 

The finite size method for j 8 c splits the phase space of the complex Polyakov loop into an outer 
region (Deconfinement phase) and an inner (Confined phase) and determines the coupling jS at 
which point the weight of the phases is 3 to 1, because of SU (3)’s Z(3) center symmetry. A proper 
defined step function of weights Wc and Wd with s{p) = (3wc — Wd)/(3wc + Wd) see Fig. 3, then has 
a zero which is determined by a linear fit. Our finite size determinations are all done with boxes 
Ns > 3Nz except for Nz = 22 at Ns = 64. For values Nz < 16 we have checked, that the weight ratio 
method results are consistent with polynomial finite size scaling extrapolations of susceptibility 
peak positions at the first order phase transition. The weight ratio method however, and already 
for Ns > 3Nz, has smaller corrections, which only are exponentially small in the spacial size. Our 
j3c data in Table 2 exert stress on semi-analytic calculations [22, 23] but are consistent with the 
Bielefeld numbers [24, 25] as well as with Berg et. al. [26]. Deviations mostly lie within the 1 a 
error margin. 

In finite temperature lattice calculations, it is wishful to present the critical temperature in 
physical units. If there are no precise measurements of physical observables like e.g. the string ten¬ 
sion we can also opt for less physical zero temperature lattice scales. We determine two scales: the 
Sommer scale [ 8 ] denoted by ro/a and the Wilson flow s/toja as recently introduced by Fiischer 
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Figure 4: Tcffif calculated from data in Table 2 at = 8,..., 18 (data points from the right to the left). The 
fo data of [10] are used to interpolate fo by the rational form of eq. 3.1. 



[9, 10] and leave the presentation to forthcoming work [6]. In case of to and for any given j8, the 
lattice configuration is “cooled” using Runge Kutta algorithm for a classical Wilson flow until a 
stopping criterion on the energy density is reached at pseudo time to- Various different regulariza¬ 
tion’s (Wilson, Wilson improved, Clover) for the energy density can be employed. Here we will 
use Luscher’s to numbers from [10] for the range 5.96 <i8<6 .59, corresponding to the range 
8 < Vf < 18. We interpolate ln(-y/to/a) similar to \n{s/a) as in eq. 3.1 by a rational form and de¬ 
termine Tc^/to in the continuum limit a^/to —)■ 0. We display in Fig. 4 the continuum extrapolation 
for Tc^/Iq based on the Luscher data. A straight line fit, the line in Fig. 4, through all six data points 
with order corrections only, yields: 

Tcs/fo = 0.2471(3)[10] xio./ = 0.94. (3.3) 

We observe very small corrections of relative magnitude 0.004 accompanied by a small statis¬ 
tical error 0.0003. Varying the fit range, and also adding corrections to the fit, we estimate an 
systematic uncertainty of 0.0010 in the observable and thus Tc-^/Jq = 0.2471(13) is a realistic find¬ 
ing. It is to be remarked, that the scaling deviations in the amplitude product turn out to be quite 
smooth and regular. Thus, under the provision of correct and precise Iq data \ our determinations 
for At <18 appear to be of similar numerical quality. Naturally the physical temperature is more 
relevant than ^o^ unless to can be converted into another physical scale at same precision levels. 

4. Conclusion 

We have estimated the critical couplings for the Wilson action and for temporal extents up 
to At = 20 with relative precision levels <0.001. For the At interval 6, ...,22 we also provide a 
rational interpolation to \n{tc/a) in eq. 3.2 and eq. 3.3, on the same footing as for the Sommer scale 
in [21], tc being the critical time extent. It remains to be seen how the determination of the diffusion 

*The values of fo/a^ in [10] are probably affected by downward finite size effects of the order few percent, as the 
authors concede. 
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constant K in eq. 2.6 responds to more elaborate models for the spectral function and depends on 
the finite quark mass and on the continuum limit. Work in this direction is in progress 

Acknowledgements: Simulations where performed on the JARA0039 and JARA0108 ac¬ 
counts of JARA-HPC in Aachen, JUDGE/JUROPA at JSC Julich, the OCuLUS Cluster in Pader- 
born, and the Bielefeld GPU cluster. 
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